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Abstract 

We introduce a method to achieve three-dimensional dynamic localization of light. We consider 
a dynamically-modulated resonator lattice that has been previously shown to exhibit an effective 
gauge potential for photons. When such an effective gauge potential varies sinusoidally in time, 
dynamic localization of light can be achieved. Moreover, while previous works on such effective 
gauge potential for photons were carried out in the regime where the rotating wave approximation 
is valid, the effect of dynamic localization persists even when the counter-rotating term is taken 
into count. 
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The effect of dynamic localization is of fundamental importance in understanding coherent 
dynamics of a charged particle in a periodic potential. When such a charged particle is in 
addition subjected to a time-harmonic external electric field, the wavefunction of the particle 
can become completely localized jl , 2], This effect has been studied in a number of systems 


131 ]. and has been demonstrated in experiments involving Bose-Einstein condensate or 
optical lattices 00 - 

Localization of photon, especially in full three dimensions, is of great practical and fun¬ 
damental importance for the control of light 12]. Hence dynamic localization of photon 
is of significance as well. Photon is a neutral particle, thus there is no naturally occur¬ 
ring time-harmonic electric field that couples to photon. To achieve dynamic localization of 
photon, one therefore needs to synthesize an effective electric field. Up to now, extensive 
experimental and theoretical works has focused on light propagation in a waveguide array, 
where the effect of dynamic localization manifests by analogy as the cancellation of diffrac¬ 


tion when the array is modulated in space along the propagation direction 18 


231. There 


has not been however any demonstration of a true three-dimensional localization of light in 
a photonic structure undergoing time-dependent modulation. 

In this letter, we show that the concept of photonic gauge potential provides a mechanism 
to achieve dy namic localization of light in full three dimensions. It has been theoretically 


proposed 


24] and experimentally demostrated [25H27] that when the refractive index of a 


photonic structure is modulated in time sinusoidal 


to an effective gauge potential for photon states 28 


H ' he phase t 
. Ref. |24 


re modulation corresponds 


31] utilized this correspon¬ 


dence to create a spatially inhomogeneous, but time-invariant gauge potential distribution, 
in order to study effects associated with an effective magnetic field for photons, including 


the photonic Aharnov-Bohm effect [24|, |25|, 127], and the photonic analogue of the integer 


quantum hall effect [28]]. In contrast, here we create a gauge potential that is spatially 
homogenous or periodic, but temporally varying. We show that such a time-dependent 
gauge potential naturally leads to a time-varying effective electric field for photons, which 
can be used to create three-dimensional dynamic localization of light. Moreover, while Ref. 


24 


28 


31] have only considered the regime where the rotating wave approximation is valid, 


here we show that such dynamic localization persist even when the counter-rotating term is 
taken into account. 

We start with the same model system as discussed in details in Refs. 


28 


29] . consisting 
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of either a one-dimensional or three-dimensional photonic resonator lattice as shown in Fig. 
CD The lattice consists of two types of resonators (A and B ) with frequencies c o A and ujb 
respectively. The Hamiltonian of the system is 

H =uj a ^ a\a n + uj b ^2 b ] n b n + V cos (^ + + b f n a m ), ( 1 ) 

n n ( mn) 

where V cos(fit + </> mn (t)) is the coupling strength between the nearest neighbor resonators, 
fi = u>a — cob- (ftmn is the phase of the couplng strength modulation. In this paper, we will 
consider the situations where such modulation phase itself is modulated in time, and refer 
to such modulation of the phase (j) mn as the phase modulation, (A (a) and (b) are the 
creation (annihilation) operators in the A and B sublattice, respectively. 

In the limit V <C fi, the rotating wave approximation (RWA) is valid. Therefore, we can 
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FIG. 1: (color online) A one-dimensional (a) and a three-dimensional (b) photonic resonator lattice 
where two kinds of resonators with frequency uj a (red dots) and ojb (blue dots). The nearest- 
neighbor coupling is dynamically modulated and the phase of the coupling constant modulation 
itself can be time-dependent with the signs being flipped between two neighboring bonds. The 
lattice is assumed infinite in all directions. 
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( 2 ) 


simplify the Hamiltonian and rewrite it in the rotating frame 


v- V 


H = ^ j (e i,l,mn(t) clc n + e i<t>mn{t) c ] n c m ) , 

(mn) 


where c m ( n ) = e ia;A ( s ha m (6 n ). In general such a system has a dynamic effective gauge held 
281 

= lmn(j>mn(t)/a , (3) 


where Z mn is a unit vector and a is the distance between two near-neighbor sites. Here 
however, we choose the modulation phases such that in Eq. (l2lh all bonds along the same 
direction have the same phase, e.g. all bonds along the ^-direction has the same phase 
4> x (t). In the three-dimensional case, 4> y (t) and 0,(f) are similarly defined. Since the phases 
are uniform in space, the system has zero effective magnetic held. 

We now show that with a proper choice of the time-dependency of these phases, we can 
achieve dynamic localization. As an illustration we consider the one-dimensional case in 
some details. The three-dimensional case naturally follows. In the one-dimensional case, as 
an intuitive analysis, we can write the Hamiltonian d2j) in the wavevector space (fc-space) 


H = Y^VcW cos l^ a - 


(4) 


Hence the system has an instantaneous photonic band structure co(k x ) = V cos [k x a—<f> x (t)\ = 
V cos[(A; x — A x )a\. The instantaneous bandstructures at three different values of 4> x are shown 
in Figure [2l The effect of a spatially uniform photonic gauge potential is a shift of the 


bandstructure in fc-space ( 29, |31]). Since the structure maintains translational invariance, 


the wavevector k x is a conserved quantity throughout the modulation process. The group 
velocity of the wave packet of the photon with wavevector k x is given by 

du{k x ) 


Vg{k x ) = 


dk x 


= —Vasm[k x a — 4> x (t)\. 


(5) 


At different values of <fi x , the group velocity at the same wavevector can have either positive 
or negative signs. (Figure [2]). 

To demonstrate dynamic localization, we choose a phase modulation of the form <f> x (t) = 
a x coswhere a x and com are the amplitude and the frequency of the phase modulation, 
respectively. Thus, the average group velocity over one phase modulation period is 


( v g {k x )) = —Va sm(k x a) J 0 (a . x ), 


( 6 ) 
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FIG. 2: (color online) The band diagram for the one-dimensional lattice as shown in Figure [TK a), 
with (a) 4> x = 0, (b) <p x = 2.40483, and (c) (j> x = —2.40483. Red dashed lines indicate the slope at 
point P. 


where Jo is the zeroth-order Bessel function, by choosing a x be to a zero of Jo, (we denote one 
of such zero as a below), the average group velocity is zero for all k x . Thus all wavepackets of 
the system become localized, signifying the presence of dynamic localization. Importantly, 
the condition for dynamic localization here is related to the strength of the phase modulation, 
and independent of the phase modulation frequency ujm- 

We confirm the intuitive analysis above based on the instantaneous bandstructure, by 
a rigorous numerical calculation of the Floquet eigenstates of Hamiltoian in Eq. (IT]) . In 


this numerical analysis, we use the Hamiltonian of Eq. (j4j), and c 


energy e at each k x , following the same procedure as in Ref. 


• fl 


irectly compute the quasi- 


32 


33| . The resulting e as 
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FIG. 3: (color online) The quasienergies as a function of a, for the Hamiltonian of Eq. (HI) . with 
4> x (t) = acos(ujMt)- Here we choose V = 0.2a jm- Each curve corresponds to a different wavevector 
k x , in the range —ir/a < k x < it/a. 


a function of phase modulation strength a, for different k x s, are plotted in Figure [3j At 
each a, the range of the values of the quasi-energy indicates the bandwidth of the quasi¬ 
energy bandstructure. The onset of the dynamic localization corresponds to the collapse of 
the bandwidth. In Figure [3j we indeed observe the collapse of bandwidth when the phase 
modulation strength approaches each of the zero’s of J 0 . 

For the study of electronic dynamic localization, the effect of a time-varying electric 
held is typically described through the use of a spatially non-uniform scalar potential, as 


described by a Hamiltonian m 
H = 


Q 

v 


Y — (<4c n + c f n c m ) - Y nauj M sin(uj M t)clc n , 


(7) 


( mn) 


In contrast, we have used a vector potential that is spatially periodic. Our Hamitonian of 
Eq. ([2]) is in fact equivalent to (Eq. (J7D) by a gauge transformation: 


I'b) = X^4|o> ->■ 1^) = = Y Vne%6n ^)i (8) 

n n n 

where |T) satishes the Schrodinger equation | H/) = H |T) or iv n = 

y_ [ e -iacos(u M t) Vn+i _)_ e *acos(w A f*) l , re _ 1 j _ With a gauge choice of 9 n = — na cos(uj M t), 
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the gauge-transformed state '!') then satisfies 

= E^ nc nl°) - E^^^l 0 ) 

n n 

= ^ E [e _ * acos( ^ Mt) u n+ i + e iacos ^ Mt ^v n -i\ e*"4|0) - E^^e^4|0) 

n n ( 9 ) 

= y E i Vn +i(? 0n+1 + Vn-xe 1 ^- 1 ] 4|0) - E naco M sin(cu m t)u n e ie "4|0) 

n n 

= m *>, 


where H is given in Eq. (0. Therefore, the two Hamiltonians of Eqs. ([2D and 0) are 
indeed equivalent to each other, as they are related by a gauge transformation. Similar 
gauge transformation has been used in the study of waveguide array [18]. Certainly, a 
time-varying gauge potential for an electron is related to an electric field applied on the 
electron. Here we have shown that a time-varying effective gauge potential for a photon also 
analogously produces an effective electric field applied on the photon. 

Unlike the waveguide array approach, where the effect of photonic dynamic localization 
manifests through an analogy as the cancellation of diffraction in a static structure, in our 
approach here one can directly achieve dynamic photon localization in all three dimensions. 
We consider the Hamiltonian of Eq. ([ID for the three-dimensional lattice as shown in Figure 
mb). We choose the phase modulation (j) Xi y, z (t) = a cos(o;m4 The intuitive derivation of 
dynamic localization condition (Eqs. (ED-(ED) can then be straightforwardly generalized to 
full three-dimension. Full three-dimensional dynamic localization is achieved provided that 
the modulation strength above is chosen to be a zero of the Jo, for all choices of phase 
modulation frequency u>m- 

Similar to the one-dimensional case, the intuitive derivations for dynamic localization 
for three-dimension can be confirmed by a rigorous Floquet analysis showing band collapse. 
Instead, here we provide the evidence of full three-dimensional dynamic localization, by a 
direct simulation of photon dynamics in a 40a x 40a x 40a three dimensional lattice. The 
simulation is done by solving the coupled-mode equation [30] 


id\V{t))/dt = H{t)\V{t)). (10) 

Here |T) = [H m n m (t)a| n + u n (f)4] |0) gives the photon state with the amplitude at site 

m{n ) described by u m ( n )(t). H{t) is the time-dependent Hamiltonian of Eq. (PQ). The initial 
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wave packet of the photon at t — 0 has the form |T(0)) = Yl v =xy z ex Pl~(v~ Vo) 2 / w2 
where (xo,yo,zo) is the center of wave packet with waist w. The results are plotted in 
Figure 0J In the absence of phase modulation, Figure HJ^a) shows the initial wave packet of 
the photon. The wave packet propagates freely in the space with time and reaches to the 
corner of the lattice at t — 1.25 a/c (see Figure 01(b)). In contrast, in the presence of phase 
modulation with a choice of the amplitude a = 2.40483 and frequency c Om = 1-5 c/a, the 
wave packet of the photon is localized near its initial position throughout the entire duration 
of the simulation. This is demonstrated in Figs. 0](c) and (d), which show the wave packet’s 
positions at t = 1.25 a/c and t = 5 a/c, respectively. The simulation here provides a direct 
visualization of the dynamic localization process in three dimensions. 

Up to this point we have used the rotating wave approximation for the Hamitonian in Eq. 
©• Previous discussions on the photonic gauge held in this Hamiltonian have all assumed 
the rotating wave approximation. On the other hand, in the experimental demonstration of 
the photonic gauge held one often uses electro-optic modulation of refractive index 26]. In 
many electro-optic modulations, the strength of the modulation, as measured in Sn/n x cuo, 
where n is the refractive index of the structure, Sn is the index change, and ujq is the 
operating frequency, can be much larger than the modulation frequency on the order of a 
few GHz, therefore it is of importance to understand the validity of gauge potential concept 
beyond the rotating wave approximation. Here we show that the dynamic localization effect 
persists even in the regime where the rotating wave approximation is not valid. 

We provide the results in one-dimension. The generalization to three-dimension is 
straightforward. For the treatment beyond the rotating wave approximation, we again 
starts by providing an intuitive treatment based on the instantaneous band-structure, we 
then confirm the intuitive treatment through an exact numerical analysis of the Floquet 
bandstructure. The Hamiltonian (0]) can be written in fc-space as: 


( 11 ) 


H = £( UAa\ x a kx + u B b\b kx ) + E' / 4>.x 

kx kx 

(e* nt cos [k x a + <f) x (t)] + e~ lQt cos [k x a — 4> x {t )]) + h.c.. 
Perform the transformation c kx = c lLJA(B)t a kx (b kx ), we obtain 

H = Vc kj c kx {cos [k x a — (j> x {t)\ + cos(2fit) cos^a + (p x (t )]} 


( 12 ) 




FIG. 4: (color online) Propagation of a photon wavepacket in a 40a x 40a x 40a three dimensional 
lattice, (a) The initial condition at t = 0, with xq = yo = zo = 20 a and k x = —k y = —k z = —1.283 
a -1 ; (b) The wave packet at t = 1.25 a/c with no phase modulation; (c) and (d) The wave packet 
at t = 1.25 a/c and t = 5 a/c, respectively, with phase modulation. The parameters of the phase 
modulation are a = 2.40483 and ujm = 1-5 c/a. The coupling strength between the resonators is 
V = 2.47T c/a. 

We notice that the first term is the same as Eq. (J3j) and the second term is the counter¬ 
rotating term. From Eq. (1T2|1 we can straightforwardly obtain the instantenous bandstruc- 
ture and hence the instanteuous group velocity at a wavevector k xi since the Hamiltonian in 
the presence of counter rotating term is still periodic in real space. Again, assuming that the 
modulation phase (j) x = a cos(u>Mt). The average group velocity over one phase modulation 
period (2iy/uj m ) is 

(v g (k x )) =—Vasm(k x a)Jo(a) — Va —— / dt cos(2Qt) sm[k x a — acos(u)Mt)}■ (13) 

2k J o 
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( v g(k x )) CR = Vax 


To facilitate analytic calculation, we assume that 

2 12 = nuiMi ( 14 ) 

where n is a positive integer, the second term in Eq. (fT3]L denoted as (v g (k x )) cr, can be 
calculated analytically as 

_l)m+ i J n («) n = 2m 

( 15 ) 

(—l) m cos(k x a) J n (a) n = 2m + 1 

By choosing a = 2.40483, which corresponds to Jo (a) = 0, the first term in Eq. (1131) van¬ 
ishes. And the correction due to the second term can be made arbitrarily small by choosing 
a sufficiently large n in Eq. i.e. by choosing the phase modulation frequency to be suffi¬ 
ciently small as compared to the frequency of coupling strength modulation. Thus, dynamic 
localization can still be accomplished in the regime where rotating wave approximation no 
longer applies. This result can be straightforwardly generalized to there-dimension. Three- 
dimensional dynamic localization should occur, when 212 = tiujm, provided that all bonds 
along each direction has the same phase 4> x . ytZ (t) = acos(cjM^) with the phase modulation 
amplitude a being a zero’s of J 0 . 

We confirm the intuitive analysis above by calculating the Floquet bandstructure in the 
case where V = 0.20, and hence the rotating wave approximation is no longer valid (blue 
lines in Figure [5]), and by comparing such calculations to the prediction of the range of 
quasi-energies with rotating wave approximation. (Red lines in Figure [5]). Figure [5](a) 
shows the case with 12 = 2Introducing the counter rotating term indeed modifies the 
bandstructure. Nevertheless, the bandwidth still collapses near a phase modulation strength 
of a = 2.40483. Thus, dynamic localization still occurs in this system beyond the rotating 
wave approximation. Figure [5](b) shows the case with 12 = 4 ujm- Comparing Figures G3(a) 
and[5](b), we observe that the discrepancy in the bandstructures between the cases with 
or without rotating wave approximation becomes smaller as ujm is reduced, in spite of the 
fact that with V = 0.212 we are significantly outside the regime where the rotating wave 
approximation is valid. This observation is consistent with the analytic results derived above 
based on instantaneous bandstructure. 

Experimentally, effective gauge held for photons have already been experimentally ob¬ 
served using two modulators [26]. The demonstration of the theoretical proposal here re¬ 
quires further integration of larger numbers of modulators. The experimental feasibility of 
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FIG. 5: (color online) Quasi-energies as a function of phase-modulation strength a, for the Hamil¬ 
tonian of Eq. m, with 4> x (t) = acos(u)Mt)- V = 0.2fi. (a) II = 2u jm- (b) = 4u jm- Each blue 
curve corresponds to a different wavevector k x . in the range of —n/a < k x < ir/a. Dashed red line 
is the envelope for the same Hamiltonian, but calculated using the rotating wave approximation. 


such integration has been discussed in Ref. 


28j. While for illustration purpose we have 


focused on a photonic gauge potential through the use of temporal refractive index modu¬ 
lation, the concept here should be relevant for other proposals of photonic gauge potential 
as well, include those based on magneto-optical effects (34|, |35|], as well as spin-dependent 


photonic gauge potential 


3-401 


and optomechnanical gauge potential 4jJ, [42]. In summary, 


we have shown that three-dimensional dynamic localization of light can be achieved with an 
effective gauge potential for photons. The results provide additional evidence of the exciting 
prospects of photonic gauge potential for the control of light propagation. 
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